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Modelling  of  PEM  fuel  cells  has  long  been  an  active  research  area  to  improve  understanding  of  cell 
and  stack  operation,  facilitate  design  improvements  and  support  simulation  studies.  The  prediction  of 
activation  polarization  in  most  PEM  models  has  concentrated  on  the  cathode  losses  since  anode  losses  are 
commonly  much  smaller  and  tend  to  be  ignored.  Further  development  of  the  anode  activation  polarization 
term  is  being  undertaken  to  broaden  the  application  and  usefulness  of  PEM  models  in  general. 

Published  work  on  the  kinetics  of  the  hydrogen  oxidation  reaction  (HOR)  using  Pt(h  /<  /)  electrodes 
in  dilute  H2S04  has  been  recently  reassessed  and  published.  Correlations  for  diffusion-free  exchange 
current  densities  were  developed  and  empirical  predictive  equations  for  the  anode  activation  polarization 
were  proposed  for  the  experimental  conditions  of  the  previously  published  work:  Pt(l  0  0),  Pt(l  1  0)  and 
Pt( 111)  electrodes,  Ph2  of  1  atm,  and  temperatures  of  1,  30  and  60  °C.  It  was  concluded  that  the  HOR  on 
Pt(l  1  0)  electrodes  followed  a  Tafel-Volmer  reaction  sequence. 

The  aim  of  the  present  paper  is  to  generalize  these  Tafel-Volmer  correlations,  apply  them  to  published 
data  for  Pt(l  1  0)  electrodes  and  further  develop  the  modelling  of  anode  activation  polarization  over  the 
range  of  operating  conditions  found  in  PEMFC  operation. 

Crown  Copyright  ©  2010  Published  by  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  background  for  the  present  paper  has  been  recently 
published  [1-3].  Our  long-term  aim  is  to  further  develop  our  gen¬ 
eralized  steady-state  electrochemical  model,  the  GSSEM,  which 
consists  of  a  simple,  1 -dimensional,  model  of  a  PEM  fuel  cell 
(PEMFC)  based  on 

V  =  E  +  TJ  act,  a  +  Pact,c  +  Pohmic  +  Pconc,a  +  Pconc,c  (1 ) 

The  specific  goal  of  the  present  paper  is  to  further  evaluate  pre¬ 
viously  published  work  which  studied  the  hydrogen  oxidation 
reaction  on  Pt(h/d)  electrodes  in  dilute  H2S04  and  to  develop  the 
results  into  a  modelling  approach  for  rjact,a,  the  anode  activation 
polarization,  in  PEMFCs.  The  quantitative  prediction  of  this  partic¬ 
ular  loss  is  generally  not  given  much  attention  in  a  PEMFC  model 
since  it  is,  normally,  much  less  significant  than  the  corresponding 
loss  at  the  cathode.  There  are  situations,  however,  where  the  anode 
activation  polarizations  could  become  significant  and  a  modelling 
capability  is  therefore  desirable. 


*  Corresponding  author  at:  Department  of  Chemistry  and  Chemical  Engineering, 
Royal  Military  College  of  Canada,  1 1  General  Crerar  Crescent,  Kingston,  Ontario,  K7K 
7B4  Canada.  Tel.:  +1  613  541  6000x6981;  fax:  +1  613  542  9489. 

E-mail  address:  thurgood-c@rmc.ca  (C.P.  Thurgood). 


Our  present  interest  is  the  anode  activation  polarization  related 
to  the  ‘hydrogen  oxidation  reaction’  (HOR)  for  which  the  overall 
half-cell  reaction  is 

H2  =>  2H+  +  2e_  (2) 

The  reverse  of  reaction  (2)  is  the  ‘hydrogen  evolution  reaction’ 
(HER). 

The  anode  reaction  occurs  in  the  presence  of  a  catalyst,  typically 
Pt  if  the  anode  feedstock  is  pure  hydrogen. 

There  is  broad  support  in  the  literature  for  the  HOR  on  Pt  to  pro¬ 
ceed  via  a  Tafel-Volmer  reaction  sequence,  fully  described  recently 
[3].  There  is  not,  however,  complete  agreement  on  the  kinetic 
model  as  literature  support  exists  for  both  ‘rds  Tafel-fast  Volmer’ 
kinetics  and  for  ‘fast  Tafel-rds  Volmer’  kinetics.  These  were  respec¬ 
tively  referred  to  as  ‘Mechanism  I’  and  ‘Mechanism  IT  in  our  recent 
reassessment  [3]  of  a  body  of  work  published  by  what  we  chose  to 
call  the  ‘Markovic  Group’  [4-9].  This  group  of  publications  included 
extensive  results  for  rotating  disk  electrode  (RDE)  studies  of  the 
HOR  on  single-crystal  Pt(l  0  0),  Pt(l  1  0)  and  Pt(l  1 1)  electrodes  in 
dilute  sulphuric  acid  electrolyte  at  1 , 30  and  60  °C  and  ‘corrected  to 
apn2  ofl  atm’.  These  publications  [4-9]  concluded  that  the  Pt(l  1  0) 
results  were  best  explained  by  ‘rds  Tafel-fast  Volmer’  kinetics,  our 
so-called  ‘Mechanism  I’. 

The  present  paper,  in  Section  2,  will  review  the  detailed  devel¬ 
opment  of  thermodynamic  and  kinetic  equations  for  ‘Mechanism 
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Nomenclature 

b 

Tafel  slope,  2.3 RT(aF)~'  (Vdec-1) 

B-V 

Butler-Volmer 

c 

concentration  (mol cm-3) 

Ch2 

concentration  of  dissolved  hydrogen  at  the  reaction 
interface  (molH2  cm-3) 

Cm 

concentration  of  active  Pt  sites  [in  sites  cm-2  (order 
1015)  or  moles  of  sites  cm-2  (order  10-8)] 

E 

equilibrium  EMF  of  a  cell  (V)  or  total  anode  polar¬ 
ization,  7/act,a  +  J)conc,a.  in  Refs.  [4-9]  (V) 

MOflo) 

functions  of  0  and  00  in  the  general  Butler-Volmer 
equation  (see  following  Eq.  (25)) 

F 

Faraday’s  Constant  (clbs  equ’t-1 ) 

GSSEM 

generalized  steady  state  electrochemical  model 

H 

Henry’s  Law  ‘constant’  (atm  cm3  mol-1 ) 

HER 

hydrogen  evolution  reaction  (reverse  of  Eq.  (2)) 

HOR 

hydrogen  oxidation  reaction  (Eq.  (2)) 

i 

current  density  (A cm-2) 

io 

exchange  current  density  (A cm-2) 

k 

chemical  rate  constant 

^ads 

rate  constant  for  the  2-site  dissociative  chemisorp¬ 
tion  of  a  hydrogen  molecule  [cm5  s-1  (mol  of  vacant 

Pt  sites)-2] 

^des 

rate  constant  for  the  hydrogen  molecule  desorp¬ 
tion  reaction  [(molH2)cm2  s-1  (mol  of  occupied  Pt 
sites)-2] 

^et.fwd 

forward  rate  constant  for  the  electron  transfer 
(Vogel)  reaction 

ket.rev 

reverse  rate  constant  for  the  electron  transfer 
(Vogel)  reaction 

Kads 

adsorption  equilibrium  constant  for  chemisorbed 
hydrogen  molecules  (also  equal  to  fcads/fc des) 
[cm3  (mol  H2  )-1 ,  the  reciprocal  of  the  units  of  Ch2  ] 

^etCHso+ 

0  equilibrium  ‘constant’,  representing 

{/<_etCH3o+  0(/<et)— 1 },  for  the  electron  transfer 
reactions,  Eqs.  (8)  and  (9) 

Mechanism  I  the  ‘rds  Tafel-fast  Volmer’  HOR  process 

n 

number  of  electrons  being  transferred  for  one  act  of 
the  overall  reaction 

na 

number  of  electrons  being  transferred  ‘after’  the  rds 

nb 

number  of  electrons  being  transferred  ‘before’  the 
rds 

P 

partial  pressure  of  a  gas  component  (atm) 

r 

chemical  reaction  rate 

rds 

rate-determining-step  in  the  reaction  sequence  at 
the  anode 

RDE 

rotating  disk  electrode 

V 

cell  voltage  (V) 

Greek  letters 

a 

transfer  coefficient 

P 

symmetry  factor 

p 

polarization  (i.e.  overvoltage  or  loss)  (V) 

0 

fractional  surface  coverage  (generally  of 
chemisorbed  hydrogen  atoms) 

V 

stoichiometric  coefficient  in  Eq.  (5)  (number  of 
times  that  the  rds  must  take  place  for  the  overall 
reaction  to  occur  once) 

Subscripts 

a 

anode  or  activity 

act 

activation 

ads 

adsorption 

b 

bulk 

c 

cathode 

cone 

concentration  (relating  to  mass  transfer  losses) 

des 

desorption 

et 

electron  transfer  (as  in  the  Volmer  reaction) 

expl 

experimental 

fwd 

forward 

h2 

hydrogen  molecule 

H 

hydrogen  atom 

H30+ 

hydrated  proton  (i.e.  H+:H20) 

i 

interface 

I 

Mechanism  I 

MG 

Markovic  Group  (authors  of  Refs.  [4-8]) 

ohmic 

relating  to  ohmic  (i.e.  iR)  losses 

0 

at  zero  polarization  and  zero  net  current  condition 
(i.e.  at  equilibrium) 

pred 

predicted 

rev 

reverse 

sat 

saturated 

T 

Tafel  reaction 

V 

Volmer  reaction 

vac 

vacant 

I’,  the  ‘rds  Tafel-fast  Volmer’  HOR  process.  Then,  in  Section  3, 
these  equations  will  be  applied  to  the  further  analysis  and  correla¬ 
tion/recorrelation  of  some  published  results  for  the  HOR  on  Pt(l  1  0) 
electrodes  and  evaluation  of  the  Tafel-Volmer,  T-V,  parameters. 
Section  4  will  then  deal  briefly  with  the  generalization  of  these 
correlations  to  values  of  pu2  and  temperatures  that  are  more  rep¬ 
resentative  of  PEMFC  modelling  applications. 

2.  Development  of  mechanistic  equations  for  a 
‘Tafel-Volmer’  HOR  reaction 

2.1.  Introduction 

The  basic  theories  of  activation  polarizations  for  the  HOR  and, 
especially,  overvoltages  for  the  HER  have  been  established  for 
decades.  This  section  will  summarize  the  development  of  a  correla¬ 
tion  and  modelling  approach  for  the  anode  activation  polarization 
that  is  as  mechanistic  as  possible  yet  as  simple  to  understand  and 
apply  as  possible.  The  goal  is  to  develop  a  relationship  linking  the 
anode  activation  polarization  to  the  current  density  and  the  other 
operating  parameters  of  the  fuel  cell. 

The  literature  dealing  with  both  the  HOR  and  the  HER  is  vast 
but  there  appears  to  be  general  agreement  on  the  main  features 
of  the  process.  Much  of  the  following  development  has  benefited 
from  the  previous  work  of  Parsons  [10],  Austin  [11,12],  Bockris 
and  Reddy  [13],  Gileadi  et  al.  [14],  the  ‘Markovic  Group’  [4-9], 
Springer  et  al.  [15],  Camara  et  al.  [16],  Breiter  [17]  and  Lasia 
[18]. 

2.2.  Fundamental  steps  in  the  HOR 

Various  possible  steps  in  the  overall  HOR  have  been  suggested, 
one  combination  being  the  so-called  Tafel-Volmer  sequence’. 
Using  ‘M’  to  represent  an  active  Pt  atom  and  ‘e-’  to  represent  an 
electron,  the  steps  in  this  sequence  are: 

(i)  H2  +  M^M.H2 

(ii)  M.H2  +  M  2M.H 

with  (i)  +  (ii),  the  ‘Tafel  step’,  often  simply 
(hi)  H2  +  2M^Had  +  Had 
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(iv)  M.H  +  H20  H30+  +  M  +  e" 

(v)  h3o+^h+  +  h2o 

with  (iv)  +  (v),  the  ‘Volmer  step’,  often  simply 

(vi)  Had^H+  +  e- 

2.3.  The  physical  picture 

At  the  gas/electrolyte  interface,  a  hydrogen  gas  partial  pressure, 
Ph2  *  is  in  equilibrium  with  concentration  cH2,Sat  of  dissolved  hydro¬ 
gen  according  to  a  Henry’s  Law  equilibrium  expression.  Dissolved 
hydrogen  will  undergo  bulk  diffusion  through  the  electrolyte  to  the 
vicinity  of  the  anode,  delivering  a  bulk  concentration  cH2,b  at  the 
electrolyte/reaction  boundary  layer  interface.  If  the  mass  transfer 
through  the  electrolyte  is  very  rapid,  or  if  it  has  been  accounted  for 
by  subtracting  the  anode  concentration  polarization  from  the  total 
anode  polarization,  cH2,b  can  be  taken  as  equal  to  cH2,Sat-  The  dis¬ 
solved  hydrogen  concentration  at  the  reaction  interface,  Ch2j  (or, 
simply,  Ch2  )  may  be  different  than  cH2?b-  Likewise,  the  proton  con¬ 
centration  at  the  reaction  interface  will  be  denoted  as  cH3q+  and 
that  at  the  interface  between  the  electrical  boundary  layer  and  the 
bulk  electrolyte  will  be  cH30+  b-  Again,  cH3q+  may  have  a  value  dif¬ 
ferent  than  cH3q+  b-  Just  as  cH2?b  is  a  function  of  factors  external  to 
the  reaction  boundary  layer,  so  will  cHb0+  b  depend  on  factors  exter¬ 
nal  to  the  electrical  boundary  layer,  e.g.  the  diffusion  rate  of  the 
hydrated  protons  through  the  electrolyte.  At  the  zero-current/zero- 
polarization  situation,  cH2  and  cH30+  at  the  reaction  interface  will 
be  denoted  as  Ch2,0  and  cH3q+  b,  respectively. 

2.4.  Kinetic  expressions  for  the  Tafel  and  Volmer  reactions 

Kinetic  expressions,  initially  in  terms  of  the  chemical  reaction 
rate  ‘r’,  for  the  above  reactions  can  initially  be  written  in  their  sim¬ 
plest  form  as  follows,  essentially  using  the  nomenclature  in  Breiter 
[17]  and  applying  Langmuir-Hinshelwood  methodology: 

For  the  Tafel,  ‘adsorption/desorption’,  reaction  involving  1 
molecule  of  H2: 

rT  —  rT,ads  —  *T,des  —  ^T,fwd^H2(l  —  ^h)  —  ^T,rev^n  (3) 

For  the  Volmer,  ‘electron  transfer  or  ionization’,  reaction,  which 
would  have  to  occur  twice  (i.e.  a  stoichiometric  number  of  2)  per 
molecule  of  H2: 

=  rVj fwd  -  *v,rev  =  /<v,fwd^H  exp[av,fWd^(^T)_1  ] 

— /<V,rev(^  —  ^h)Uh+ exp[— Qfv,revPT(i?T)  ]  (4) 

In  Eqs.  (3)  and  (4),  the  rate  constants,  /q,  will  include  some 
parameter  to  represent  the  catalyst  concentration. 

The  transfer  coefficient,  a,  has  recently  been  discussed  [2,3]  and 
can  generally  be  expressed  as  a  function  of  /3,  the  symmetry  factor. 
Based  on  earlier  work,  the  following  expression  was  proposed  by 
Gileadi  et  al.  [14]  with  the  symbols  and  subscripts  as  defined  in  the 
Nomenclature: 

u  =  [jS(n-nb-na)  +  nb](v)“1  (5) 

In  kinetic  expressions  involving  catalyst  concentration,  the  fraction 
of  the  Pt  sites  that  are  vacant  is  often  required.  This  leads  to  the 
introduction  of  the  common  Langmuir-Hinshelwood  parameter, 
6V ac,  the  fraction  of  the  active  sites  that  is  vacant  and  therefore  able 
to  participate  in  the  chemisorption/dissociation  reactions  on  each 
electrode.  For  a  feed  that  is  free  of  impurities,  this  should  simply 
depend  on  the  fraction  of  the  sites  occupied  by  adsorbed  H  atoms, 
^Hads*  i-e-  ^vac  =  1  -  ^Hads  •  Since  the  present  paper  considers  that  ‘H’ 
(i.e.  Hads)  is  the  only  adsorbed  species,  0  will  be  understood  to  mean 

0Hads- 


Mechanistic  equations  for  the  general  Tafel-Volmer  reaction 
sequence  will  now  be  developed  in  detail,  enlarging  on  and  modify¬ 
ing  the  Eqs.  (3 )  and  (4)  proposed  by  Breiter  and  utilizing  the  detailed 
treatment  in  Section  5.9.3  of  Austen  [11].  The  special  case,  Mecha¬ 
nism  I,  ‘rds  Tafel  followed  by  fast  Volmer’,  will  then  be  developed 
from  these  general  equations. 

2.5.  Tafel-Volmer  equations  for  the  HOR 

2.5 A.  Introduction 

The  following  development  of  equations  to  describe  the  HOR 
is  based  on  chemisorption  according  to  a  Langmuir  isotherm,  i.e. 
active  Pt  sites  energetically  homogeneous  and  no  interactions 
between  adsorbed  species.  Although  ‘ideal’  Langmuir  adsorption 
is  often  assumed  because  of  the  relative  simplicity  of  the  result¬ 
ing  kinetic  equations,  the  Markovic  Group,  in  fact,  did  study  the 
chemisorption  of  hydrogen  on  the  various  Pt(h/d)  crystal  faces. 
Their  study  [5],  in  conjunction  with  the  earlier  work  of  Vetter  [19], 
concluded  that,  although  the  Pt(l  1  l)data  were  best  fitted  by  anon- 
ideal  Frumkin  isotherm,  “the  HOR  on  Pt(l  1  0)  at  low  overpotentials 
followed  application  of  the  Langmuir  ideal  adsorption  isotherm  of 
Hads  and  the  ideal  dual-site  form  of  the  Tafel-Volmer  sequence, 
the  atom-atom  recombination  step  (i.e.  the  reverse  Tafel  reaction) 
being  the  rds.”  In  addition,  they  concluded  for  Pt(l  1  0)  “that  the 
kinetics  and  mechanisms  of  the  HER  and  HOR  should  be  the  same”. 

The  ‘nF  conversion  from  a  chemical  rate  equation  (mole- 
based)  to  an  electrochemical  rate  equation  (current-based)  is  now 
added.  Parameter  CM  represents  the  active  catalyst  concentration 
at  the  reaction  interface.  Typically,  for  low-pressure  and  near  ideal 
behaviour,  activities  can  be  replaced  by  the  corresponding  con¬ 
centrations.  The  <3h2  in  Eq.  (3)  will,  therefore,  be  replaced  by  the 
hydrogen  partial  pressure,  pH2 ,  or,  more  appropriately  for  a  liquid- 
phase  reaction,  cH2,  the  concentration  of  dissolved  hydrogen  in  the 
electrolyte  at  the  reaction  interface.  Similarly,  the  aH+  in  Eq.  (4)  will 
be  replaced  by  cHs0+,  the  concentration  of  the  hydrated  protons  at 
the  reaction  interface. 

2.5.2.  The  Tafel  reaction  for  dissociative  chemisorption 

Eq.  (3),  with  ‘n  =  2’  and  with  the  refinements  described  above, 
leads  to  the  rates  of  adsorption  and  desorption  (per  unit  catalyst 
surface  area),  expressed  as  current  densities,  as  follows: 

il.ads  =  2FkadscH2[CM(l  -  6)]2  =  2F(kadscH2C2  )(1  -  ef  (6) 

iT,des  =  2 Fkdes(CM0)2  =  2F(kdesc2  )02  (7) 

As  suggested  in  Section  2.3,  the  concentration  of  dissolved  gas  at 
the  reaction  interface,  Ch2  ,  could  be  considered  as  an  unknown  and 
cannot,  therefore  be  evaluated  independently  from  kads  in  Eq.  (6).  It 
can,  however,  be  replaced  by  the  product  {CH2(cH2,b)_1KH2,b  with, 
as  previously  suggested,  cH2,b  ideally  equal  to  CH2,Sat-  The  empirical 
ratio  {CH2(cH2,b)_1}  could  be  a  function  of  current  density  and,  for 
simplicity,  will  be  represented  by  the  symbol  yw2 .  The  zero-current 
value  of  the  ratio,  Fh2,o>  representing  the  ratio  CH2,0(cH25b)_1  *  should 
be  unity.  If  mass  transfer  of  dissolved  H2  to  the  interface  is  suffi¬ 
ciently  rapid,  yH2  at  all  current  densities  should  be  unity. 

2.5.3.  The  Volmer  equations  for  electron  transfer 

Eq.  (4),  with  ‘n  =  V  but  with  a  stoichiometric  number  of  2, 
and  with  the  refinements  described  above,  leads  to  the  following 
expressions  for  the  current  densities  of  the  electron  transfer  (i.e. 
Volmer)  reactions: 

*'v,fwd  —  2E/<et,fwd(Qvl$)  CXp[ctfv,fwdTp(RT)  ]  (8) 

*V,rev  —  2F/<et,revCj33o+  {QviO  —  ^)l  CXp[—  dyvevFq{RT)  ]  (9) 
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The  anode  activation  polarization  is  given  by  q  although  the  sign  in 
the  exponential  terms  of  Eqs.  (8)  and  (9)  can  be  changed  depending 
on  the  desired  sign  of  the  calculated  q. 

As  for  Ch2  the  Tafel  reaction,  the  concentration  of  hydrated 
protons  at  the  reaction  interface,  cH3q+,  is  unknown  and  can¬ 
not,  therefore,  be  evaluated  independently  from  ke t.rev  in  Eq. 
(9).  A  combined  unknown,  {/<et,revCH3o+,o)  *s  therefore  defined 

and  an  empirical  ratio,  Ch3o+(ch3o+,o)_1’  is  proposed  such 
that  the  group  7<et,revCH3o+’  can  be  replaced  by  the  group 

{fcet,revCH30+0){cH30+(cH30+  0)-1}.  For  simplicity,  the  dimen¬ 
sionless  ratio  cH30+(cH30+  0)_1  will  be  represented  by  the  symbol 
Kh3o+-  At  zero  current,  cH’30+  equals  cH30+  0  so  that  Kh3o+,o  then 
becomes  equal  to  unity. 

2.5.4.  The  Volmer  reaction  if  Mechanism  I  applies 

If  the  Volmer  reactions  are  very  fast  compared  to  the  Tafel 
reactions,  the  Tafel  reactions  become  the  rate-determining  step, 
the  forward  and  reverse  Volmer  currents  are  each  very  large 
compared  to  the  Tafel  currents  and  the  Volmer  reactions  are  at 
near  equilibrium.  In  other  words,  equating  Eqs.  (8)  and  (9), 

2F/<et,fwdCM  exp[av,fwdfr?(fi'T)  1  ]0  ~  2F{/<et,revCH30+,o}CM>T[30 

+  exp[-av,revF»)(Rrr1](l-0)  (10) 

With  the  assumptions  [3]  that  aty.fwd  equals  fiy  and  c^v, rev  equals 
(1  -  fiv),  Eq.  (10)  leads  to  the  following  two  expressions: 

0  ^  {/<et,revCH3o+?0}(^et,fwd)  1  Fh30+  exP(  — Frj(RT)  1 )} 

X  [1  +  {/<et,revCH3o+)0}(^et)  Yh30+  +  exPl  — Fq(RT)  }]  (H) 

Fh30+  ~  °0  -  ^r1fcet,fwdffcet,revCH30+  or1  exp{ F^RF)-1}  (12) 

Letting  /<et,rev(ket,fwd)_1  be  represented  by  Ke t,  the  equilibrium 
constant  for  the  electron  transfer  reaction,  Eq.  (12)  becomes 

Kh3o+  «  00  ~  0)“’  exp{F»7(Rr)-1){KetCH30+  (13) 

2.5.5.  The  Tafel  reaction  if  Mechanism  I  applies 

If  the  Volmer  reaction  is  very  fast  relative  to  the  Tafel  reac¬ 
tion,  the  forward  and  reverse  Volmer  currents  will  essentially 
determine  the  value  of  0.  The  ‘slow’  Tafel  reaction  will  then  be 
the  rate-determining-step  and,  even  though  the  electrons  are  all 
produced  by  the  Volmer  reaction  (as  seen  in  Section  2.2),  it  will 
be  the  basis  for  predicting  the  net  (i.e.  observed)  current  density. 
Therefore,  from  Eqs.  (6)  and  (7): 

i  =  2F(kadscH2C2  )(1  -  6)2  -  2 F(kdesc2  )02  (14) 

2.5.6.  The  ‘zero  current,  zero  polarization ’  expressions 

The  exchange  current  density,  essentially  the  rate  of  reaction 
in  chemical  kinetics  terminology,  reflects  the  kinetics  of  the  Tafel 
reaction  at  the  60  condition.  Here  there  is  no  net  current,  the 
exchange  current  density  i0  is  equal  to  both  iads  0  and  ides  0  and 
the  equilibrium  value  of  6 ,  00 ,  is  consistent  with  the  adsorp¬ 
tion/desorption  equilibrium  given  by  the  equality  of  Eqs.  (6)  and 
(7): 

*0  —  *ads,0  —  2F(kads^H2,b^M^H2,o(^  —  ^o)  (15) 

*0  —  *'des,0  —  2F(/<desC]^)0Q  (lb) 

By  definition,  yu2,o>  representing  CH2,o(cH2,b)_\  is  unity  at  zero 
current  so  that  this  parameter  should  disappear  from  Eq.  (15). 
With  kads(^des)-1  being  replaced  by  the  adsorption  equilibrium 
constant  JCads,  the  following  equilibrium  expressions  result: 

^(l-tfo)-1  =  (iWH2,b)1/2  01) 


Oo  =  (KadsCH2,b)1/2fl  +  (KadsCH2,b)1/2f 1  OS) 

Also,  from  Eqs.  (8),  (9)  and  (1 1 ),  6q  can  be  given  by: 

Oq  ~  {KetCj^3o+  0i[l  +  i^etcH30+,oi]  (19) 

This  represents  the  maximum  value  of  6  since,  once  current  starts 
to  flow,  the  discharge  reaction  reduces  the  coverage  of  H  atoms  on 
the  surface  and  6  decreases  from  60. 

As  the  polarization  approaches  zero,  0  in  the  electron  trans¬ 
fer  expressions,  Eq.  (13),  with  yH3o+  now  also  approaching  unity, 
must  approach  the  equilibrium,  ‘zero-current’,  value,  60,  given  by 
Eq.  (17).  This  leads  to  the  requirement  that 

OoO-  OoT1  *KetcHs0+  o  (20) 

Eqs.  (17)  and  (20),  therefore,  imply  that 

(21) 

The  equilibrium  constants  /Cads  and  I<e t  should  be  functions  of  tem¬ 
perature  only  so  that,  as  indicated  by  Eq.  (21 ),  the  zero-polarization 
value  of  cH3q+  ,  cH30+?0*  should  adjust  until  it  is  approximately  equal 

to  the  ratio  t)-1  at  the  reaction  interface.  The  degree  of 

approximation  depends  on  the  relative  values  of  the  Volmer  and 
Tafel  reactions.  If,  for  example,  the  forward  and  reverse  Volmer 
reactions  are,  respectively,  infinitely  faster  than  the  forward  and 
reverse  Tafel  reactions,  the  Tafel  reactions  have  an  insignificant 
influence  on  the  value  of  0  and  the  equalities  in  Eqs.  (19)-(21) 
become  exact. 

2.5.7.  Conversion  to  the  form  of  a  Butler-Volmer  equation 
Dividing  Eq.  (14)  by  Eqs.  (15)  and  (16)  leads  to 

i(»o)_1  =  YH2(yH2,o)-'[0  -0)(1  -Oo)-']2  -  moT1]2  (22) 

where  yH2(FH2,o)_1  is  actually  simply  cH2(cH2,0r1- 

Combining  Eq.  (11)  with  Eq.  (22)  leads  to  the  following  expres¬ 
sion  linking  current  density,  i,  and  activation  polarization,  q: 

i  =  ioj[cH2(cH2, o)-1(l  -  0)2{\  -  Oo)  2 ]{1  -  {Kh3o+}2 

X  [expl-F^RF)-1)]2)  (23) 

This  is  of  the  general  form  of  a  Butler-Volmer  equation.  If  Eq.  (23)  is 
rearranged  into  the  form  of  a  Tafel  equation,  the  following  results: 

f?  =  -RF(2F)-1  ln[{yH30+r2 

X  [1  -{J(io.Tr1[CM2(cn2,or1(l -0)2O  -flo)-2]”1)]  (24) 

Austen  [1 1  ],  via  a  more  detailed  route,  derived  the  same  equations 
(in  Section  5.9.3,  Case  2A,  around  Eqs.  (5.68)-(5.71))  and  Breiter 
[17]  presented  a  similar  equation.  For  brevity,  the  function  can  be 
represented  by  the  symbol  f  and  the  expression  can  be  put  in  the 
traditional,  ‘log  base  10’,  form  so  that  Eq.  (24)  can  be  written 

r, I  =  -2.303R7X2F)-1  log[(yH30+r2)l  -  idoj/l)-1)]  (25) 

Eq.  (25)  provides  the  basis  for  the  Tafel  plot  although,  if  AH3q+  and 
fi  are  significantly  different  from  unity,  the  usual  ‘q  vs.  log  i'  plot 
will  not  give  a  very  accurate  value  of  i0> T. 

3.  Evaluation  of  Tafel- Volmer  parameters  from  published 
polarization  data 

3.1.  Introduction 

The  comparatively  few  published  results  on  the  HOR  for  Pt(  1 10) 
were  reviewed.  The  various  thermodynamic  and  kinetic  equa¬ 
tions  that  have  been  derived  above  for  the  Tafel-Volmer  reaction 
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Table  1 

Exchange  current  densities,  i0,  in  Acrrr2  for  the  HOR  on  Pt(l  1  0)  from  the  RDE 
polarization  data  of  the  Markovic  Group. 


Temperature  (K) 

lo.Note  l  X  1 03 

lo,Note2  X  103 

lo,Note3  X  103 

274 

0.65 

1.0 

0.93 

303 

0.98 

1.35 

1.45 

333 

1.35 

1.50 

1.55 

Note  1 :  As  originally  recommended  by  the  Markovic  Group  [4-9]. 

Note  2:  As  subsequently  proposed  by  Mann  et  al.  [3]  based  on  the  proposed  r](i) 
database  in  Ref.  [3]. 

Note  3:  Now  proposed  using  the  Ref.  [3]  ‘limit  method’  but  using  the  expanded  rj(i) 
database  currently  proposed  in  Table  2. 

sequence  were  applied  and,  where  possible,  values  of  the  various 
Tafel-Volmer  parameters  were  determined.  The  goal  was  the  eval¬ 
uation  of  parameter  values  for  use  in  mechanistic  equations  for  the 
prediction  of  p  (i.e.  ijact,a)  as  a  function  of  current  density  (or  anode 
activation  polarization)  and  cell  operating  parameters  for  the  HOR 
on  Pt(l  10). 

3.2.  The  original  published  results  of  the  Markovic  Group 

The  ‘Markovic  Group’  body  of  work  [4-9]  included  results  at 
274  K,  303 1<  and  333  K  with  the  most  comprehensive  results  being 
at  the  lowest  temperature  studied,  274  K.  They  reported  the  i0 
results  for  their  Pt(l  1  0)  single-crystal  RDE  measurements  that  are 
listed  in  Table  1 .  For  the  reasons  recently  presented  [3  ],  their  results 
were  reassessed  and  are  further  discussed  in  Section  3.3. 

The  catalyst  is  considered  as  a  ‘reactant’  and  appears  in  kinetic 
expressions  as  a  ‘concentration’,  symbol  Cm.  The  ‘Markovic  Group’ 
concluded  that,  in  single-crystal  studies  [8],  well-ordered  Pt(h/d) 
surfaces  could  be  produced,  each  Pt  configuration  having  a  par¬ 
ticular  Pt  atom  surface  density.  Of  present  interest  is  the  surface 
density  of  Pt(l  1  0),  which  they  reported  as  0.94  x  101 5  atoms  cm-2. 
This  gave  a  Cm,i  i  o  value  of  1 .56  x  1 0-9  mol  cm-2,  the  value  used  in 
our  subsequent  data  analyses. 

3.3.  Our  reassessment  of  the  published  data  from  the  ‘Markovic 
Group’ 

3.3.1.  Our  preliminary  reassessment  of  the  ‘Markovic  Group’  data 
Our  original  reanalysis  [3]  of  the  published  results  of  the 
Markovic  Group  began  with  the  estimation  of  mass-transfer- 
corrected  ^act(0  data  sets,  shown  in  Table  2  for  their  Pt(l  1  0) 
electrodes,  from  their  previously  published  E(i)  data  sets. 
These  ‘extracted’  rjactii)  data  sets  were  then  recorrelated  using 
Tafel-Volmer  equations  summarized  in  Section  2. 

For  Pt(l  1  0),  assuming  the  ‘Mechanism  I’  special  case  of  the 
Tafel-Volmer  equations  at  all  temperatures,  the  i0  results  listed  in 


Table  2 

PactM)  data  sets  for  Pt(110)  estimated  from  the  E(i)  polarization  data  of  the 
Markovic  Group  (expanded  from  Ref.  [3]). 


274  K 

303  K 

333  K 

i  x  103  A  cm-2 

x  1 03  V 

i  x  103  A  cm-2 

T]X  103V 

i  x  103  A  cm-2 

rjx  103V 

0.0655 

0.78 

0.10 

0.90 

0.10 

0.868 

0.131 

1.45 

0.20 

1.80 

0.19 

1.571 

0.263 

3.00 

0.30 

2.69 

0.20 

1.65 

0.491 

5.50 

0.362 

3.24 

0.30 

2.304 

0.719 

7.80 

0.40 

3.52 

0.40 

2.828 

1.0 

10.7 

0.50 

4.34 

0.414 

2.891 

1.1 

12.0 

0.60 

5.15 

0.50 

3.362 

1.2 

13.5 

0.70 

5.93 

0.60 

3.863 

1.3 

15.2 

0.80 

6.66 

0.647 

4.053 

1.4 

17.0 

0.90 

7.34 

0.70 

4.26 

1.00 

7.95 

0.80 

4.65 

0.90 

4.94 

Several  datasets  have  been  added  or  slightly  altered  in  value  since  Ref.  [3]. 


Table  1  were  recommended  from  this  reanalysis.  These  i0  results 
came  from  the  following  application  of  Eq.  (23): 

iOJ  =  lim  [i{l  -  exp(-2Fr/(RT)~: )”’  ]  (26) 

i^O 

Accompanying  this  reanalysis,  empirical  expressions  were  pro¬ 
posed  to  represent  the  variation  with  current  density  of  the  group 
[Ch2(ch2l0r1(l  -  0)2O  -  OoT2},  i.e./i,  in  Eqs.  (23)-(25). 

Before  the  kinetic  and  thermodynamic  parameters  could  be 
quantified,  a  basis  for  the  prediction  of  cH2  had  to  be  established. 

3.3.2.  Hydrogen  solubility  in  the  H2S04  electrolyte  and  the 
prediction  of  cH2 

Values  of  Ch2  are  required  to  apply  Eq.  (6)  and  the  various  expres¬ 
sions  derived  from  Eq.  (6).  For  the  0.05  M  H2S04  electrolyte  used 
by  the  Markovic  Group  and  based  on  their  statement  that  ‘data 
were  reported  as  adjusted  to  a  hydrogen  gas  pressure  of  1  atm’, 
i.e.  ‘pn2  =  1\  the  following  correlations  [1]  are  considered  to  be 
applicable  for  the  prediction  of  Cn2,SaG 

Ch2,o— 45 °c  =  1.25  X  1 0-7  exp(545T-1 )  (27) 

Ch2,45-ioo°c  =  1.19  X  10"6 expt-iyor-1)  (28) 

As  discussed  in  Section  2.3,  this  Cn2,Sat  should  be  the  same  as  cH2?b  at 
the  electrolyte/reaction  boundary  layer  interface  if  diffusion  in  the 
bulk  electrolyte  has  been  properly  accounted  for.  Further,  if  there  is 
no  depletion  of  dissolved  hydrogen  in  the  reaction  boundary  layer 
adjacent  to  the  reaction  interface,  the  reaction  interface  value,  cH2, 
should  also  equal  Cn2,Sat- 

Eqs.  (27)  and  (28)  intersect  at  45  °C  so  that  a  slight  inflection 
point  might  be  expected  at  this  temperature  in  Arrhenius  plots  of 
parameters  which  vary  with  Ch2  . 

3.3.3.  Our  present  correlation  ofPt(l  1  0)  results  by  Mechanism  I 
Subsequent  to  our  preliminary  analysis  [3],  several  additional 

p(i)  datasets  were  estimated  from  the  Markovic  Group  E(i)  datasets 
(Table  2). 

The  various  T-V  models  described  in  Section  2  have  a  large  num¬ 
ber  of  adjustable  parameters  and  constraints  to  consider  during 
correlation: 

-  the  Tafel  rate  parameters  kads  and  /<des  (or,  in  some  cases,  the 
equilibrium  parameter  /Cads); 

-  the  Volmer  rate  parameters  ket,fwd  an^  ket.rev  (or,  in  some  cases, 
the  equilibrium  parameter  I<e t); 

-  the  Kh?o+(')  or  ]/h3o+('?)  variation; 

-  the  values  of  the  transfer  coefficients  aty.fwd  and  afy.rev; 

-  the  desirability  that  the  above  rate  parameters  and  equilibrium 
parameters  show  Arrhenius  behaviour; 

-  the  desirability  that  the  parameter  values  resulting  from  the  cor¬ 
relations  predict  i0  values  that  agree  with  the  ‘Note  3’  values  in 
Table  1. 

Various  combinations  of  Eqs.  (6),  (7)  and  (1 1  )-(l  6)  were  applied 
to  correlate  the  datasets  in  Table  2.  Since  Table  2  data  were  the 
result  of  a  mass-transfer  correction  to  the  original  E(i )  data,  a  value 
of  unity  was  assumed  for  the  parameter  yH2.  The  parameter  yHs  0+, 
however,  was  not  necessarily  subject  to  the  assumption  of  unity. 
Values  of  the  transfer  coefficients,  ofy.fwd  and  afy.rev.  were  taken  to 
be  the  0.5  value  previously  found  satisfactory  [3]. 

The  most  successful  correlation  technique  utilized  a  ‘general’ 
Tafel-Volmer  spreadsheet  program  so  that  kads,  /<des,  /<et,fwd*  ^et.rev 
and  yH  0+  were  all  adjustable  parameters.  To  convert  this  to  a 
Mechanism  I,  ‘rate  determined  by  Tafel  reaction  and  0  essentially 
determined  by  the  Volmer  reaction’,  values  of  /<et)fwd  and  ket.rev  were 
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Table  3 

Correlation  of  Pt(l  1  0)  ij(i)  data  from  Table  2  using  ‘Mechanism  I’  Tafel-Volmer  equations  with  yu2  set  equal  to  1  and  with  i0  consistent  with  i0,Note3  in  Table  1. 


T(  K) 

{fCtCH3o+)0}  (1) 

/<ads  X  10-16 

/<des  X  10-9 

Xh3o+ 

l'o  X  103 

/C1/2c1/2(l) 

ads  LH2  ^  J 

0o 

274 

1.3545 

1.20 

5.98 

1  to  <1 

0.93 

1.3459 

0.5753 

303 

4.800 

13.8 

4.51 

1  to  >1 

1.45 

4.800 

0.8276 

333 

14.8533 

119.5 

3.87 

1  to  <1 

1.60 

14.809 

0.9369 

(1)  As  discussed  in  the  paper,  should  be  equal  to  {^etCH3o+,o}- 


increased  until  variations  in  /<ads  and  fcdes  no  longer  affected  the  val- 
ues  of 6.  Values  of  kads,  kdes,  the  group  {ket,revCH?0+  ^(kg,^)-1 1  and, 
when  appropriate,  yH  0+  were  then  varied  until  the  best  agreement 
of  ipred  and  iexp  was  obtained  for  the  datasets  for  a  particular  tem¬ 
perature.  Typically  this  meant  an  average  ipred(*exp)-1  of  unity,  an 
i0  value  consistent  with  the  i0,Note3  values  in  Table  1,  and  the  best 
‘ipred  vs.  /eXp *  parity  plot. 

The  data  were  correlated  with  five  major  constraints  in  mind: 

(1)  The  resulting  set  of  parameter  values  should  predict  values  of 
i0  that  are  in  close  or  reasonable  agreement  with  the  i0,Note  3 
values  in  Table  1. 

(2)  The  resulting  sets  of/<ads,  /<des  and  {/<-etCH30+  0(/<et)_1}  parameter 
values  should  show  good  Arrhenius  behaviour. 

(3)  The  values  of  >'h3o+(z)  at  the  3  experimental  temperatures  can 
be  assumed  to  have  an  ideal  value  of  unity. 

(4)  The  final  values  of  {KadscH2b)^^2  and  {/CetCH30+>0}  should  be  in 
close  agreement. 

(5)  In  addition,  if  Xh3o+(1)  zs  a  physico-chemical  parameter  and  is 
not  equal  to  unity,  the  final  values  should  show  a  variation  with 
temperature  that  was  reasonable. 

Several  correlation  approaches  were  applied: 

(i)  one  that  placed  most  emphasis  on  constraints  (1)  and  (5) 
above; 

(ii)  one  that  placed  most  emphasis  on  constraints  (2),  (3)  and  (5) 
above,  and 

(iii)  one  that  placed  additional  emphasis  on  constraint  (5)  above. 

It  was  found  that  all  the  above  constraints  could  not  be  com¬ 
pletely  met.  The  best  correlation  results  are  summarized  in  Table  3 
and  plotted  in  Fig.  1.  Table  1  i0,Note3  values  were  met  at  274  K  and 
303  K  while  an  i0  value  3%  higher  than  the  i0,Note3  value  gave  the 
best  results  at  333  K. 


Anode  activation  polarization^  gct  (V) 

Fig.  1.  Calculated  values  of  do  and  6  using  proposed  ‘Mechanism  I’  equations  and 
the  Tafel-Volmer  parameter  values  in  Table  3. 


3.3.4.  Correlation  ofTafel  parameters  kads  and  kdes 

The  proposed  preliminary  Arrhenius  correlations  for  the  adsorp¬ 
tion  and  desorption  rate  constants  in  Table  3  are,  respectively: 

kads  «  2.37  x  1027  expl-y^CT)”1)  (29) 

kdes  «  5.24  x  K^expiSeOCrr1}  (30) 

These  correlations  predict  the  experimental  values  of  fcads  and  /<des 
in  Table  3  to  within  about  ±21/2%.  This  is  quite  reasonable  for  many 
kinetic  applications  but  not  really  good  enough  when  the  current 
density  at  the  anode  is  low  and  when  the  forward  and  reverse  cur¬ 
rent  densities  can  be  relatively  close  in  value.  For  this  reason,  since 
the  probable  uncertainty  in  the  parameter  values  at  the  3  experi¬ 
mental  temperatures  is  similar,  separate  correlations  are  proposed 
below  for  the  2  temperature  regions,  ‘below  303  K’  and  ‘above 


303  K’: 

kads, <303  =  1-55  x  1027  exp{-7013(7r1)  (31) 

kdes, <303  =  3.12  x  108  exp{809(r)_1)  (32) 

kads, >303  =  3.83  X  1027exp{— 7287(7T1)  (33) 

kdes, >303  =  9  06  X  108  exp{486(r)_1)  (34) 

3.3.5.  Correlation  ofTafel  parameters  I<ads  and  1  KadscH2 1'^2 

With  Kads  defined  as  the  ratio  /<ads(^des)_1*  the  following  is 
implied  by  Eqs.  (31)-(34): 

Kads, <303  =  4.97  X  1018exp{— 7822(7T1)  (35) 

K ads, >303  =  4.227  x  1018  exp{-7773(7T 1 )  (36) 


Either  by  combining  Eqs.  (35)  and  (36)  with  Eq.  (27)  or  (28),  or 
by  using  /<ads  and  /<des  values  from  Table  3,  values  of  the  equilib¬ 
rium  group  {/CadsCH2 11/2  can  be  obtained.  These  are  listed  in  Table  4. 
Correlations  for  {/‘CadsCH2 11^2  are  not  proposed  since  the  above  cor¬ 
relations  for  JCads  have,  of  necessity,  been  based  on  an  inflection  at 
303  K  while  the  published  cH2  correlations,  Eqs.  (27)  and  (28),  are 
based  on  an  inflection  at  45  °C,  i.e.  318  K. 

As  implied  by  Eq.  (21 ),  these  proposed  {/CadscH2 11/2  values  should 
also  give  reasonable  values  of  {I<etCH3o+  0}.  This  is  covered  in  the 
following  section. 

3.3.6.  Correlation  of  Volmer  parameter  {KetCH  0+  0} 

The  experimental  values  of  this  parameter  are  listed  in  Table  3. 
They  are  reasonably  fitted  by  the  following  Arrhenius  expression: 

(KetCH3o+o)  *  1.00(106)exp{-3704(r)_1)  (37) 


Table  4 

Tafel  parameters  I<ads  and  (/CadscH2  )1/2  for  Pt(l  1  0). 


T(K) 

kads(kdes)-’xl0-7(l) 

Ch2  X  107  (2) 

(kadsCH2  f/2 

274 

0.199 

9.136 

1.348 

303 

3.055 

7.552 

4.803 

333 

30.80 

7.142 

14.83 

(1)  From  Eqs.  (35)  and  (36). 

(2)  From  Eqs.  (27)  and  (28). 
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As  for  the  parameter  {I<a dscH2}1//2  in  the  previous  section,  the  45  °C 
inflection  point  can  be  assumed  giving: 

{KetCH3o+, o)0_45„c  =  746  x  105  exp{-3620(7T1}  (38) 

(KetCH30+,0}>45oc  =  2.41  x  loSexpf-sggsCT)-1)  (3g) 

The  two  sources  of  {I<etCH  0+  0)  values>  experimental  values  from 
Table  3  and  predicted  {I<etCH30+  0}  values  from  Eqs.  (38)  and  (39), 

and  the  predicted  {KadscH2}^^2  values  in  Table  4,  are  in  very  good 
agreement  as  expected  from  Eq.  (21 ). 

3.3.7.  Calculation  of  0 0  and  0 

Values  of  0o  are  estimated  from  Eq.  (19),  listed  in  Table  3,  and 
plotted  in  Fig.  1 .  As  expected,  and  consistent  with  Eqs.  (3 1  )-(36),  the 
values  of  0o  increase  with  temperature  and  reach  levels  well  above 
the  usual  ‘low-coverage’  referred  to  with  Langmuir  adsorption. 

Values  of  0  were  evaluated  via  Eq.  (11)  during  the  correlation 
process  and  representative  values  are  also  plotted  in  Fig.  1.  These 
6{p)  plots  are  of  relatively  low  slope,  possibly  indicating  that  the 
variation  in  0  is  so  low  that  non-Langmuir  behaviour  should  not 
show  up.  This  6{p)  variation  is  examined  further  in  Section  3.3.8. 

3.3.8.  Correlations  and  parity  plots  at  current  densities  above  zero 

In  general,  correlation  at  increasing  current  densities  could 
not  be  attained  if  yU3o+  was  constrained  to  a  constant  value  of 
unity.  As  listed  in  Table  3,  yH3 0+  had  a  value  at  or  near  unity 
at  low  current  density  but  generally  required  a  non-unity  value 
as  current  density  increased  if  a  good  parity  plot,  ipred  vs.  ieXp, 
was  to  be  obtained.  Graphed  values  of  Kh3o+  did  not  show  a 
consistent  variation  with  current  density,  with  r\  or  with  tem¬ 
perature  so  that  simple  correlation  was  not  possible  or,  probably, 
meaningful. 

However,  for  production  of  a  parity  plot,  extrapolation  of  the 
polarization  curves  to  higher  current  densities  or  for  prediction  of 
polarization  curves  at  other  operating  conditions,  beginning  the 
calculation  with  a  value  of  p  and  then  calculating  the  resultant  cur¬ 
rent  density  was  the  simplest  approach.  An  empirical  parameter, 
either  yH3o+  or  0* was  therefore  required  as  a  function  of  p. 

It  became  clear  that  0(p)  was  the  best  semi-empirical  function 
to  bring  in  any  departure  of  yU3 0+  from  unity,  the  following  expres¬ 
sions  representing  the  fits  of  the  0  values  required  for  the  kinetic 
parameter  values  summarized  in  Table  3.  It  was  evident  that  the 
linear  6{p)  expression  for  303 1<  was  reasonably  accurate  but  the 
linear  6(1 7)  expressions  for  274  K  and  333  K  showed  deviation  in 
opposite  directions  at  higher  values  of  current  density.  An  addi¬ 
tional  small  term  is,  therefore,  proposed  below  for  the  274  K  and 
the  333  K  empirical  expressions: 


$274  K 

«  6>o,274K  -  13.077  +  3.g3  X  10V35 

(40) 

$303  K 

=  $0,303K  -  5.457? 

(41) 

$333  K 

«  $0,333 K  -  2.377  -  1.86  X  105773  35 

(42) 

Parity  plots  using  the  experimental  results  were  then  produced  as 
follows: 


•  Eqs.  (31  )-(34)  were  applied  to  obtain  values  of  /<ads  and  /<des; 

•  Eqs.  (27)  and  (28)  were  applied  to  obtain  predicted  values  of  cH2  ; 

•  Cm  was  equal  to  1.56  x  10-9; 

•  Eqs.  (40)-(42)  were  applied  to  the  experimental  values  of  p  in 
Table  2  to  obtain  predicted  values  of  0  for  each  experimental  data 
set; 

•  Eq.  (14),  with  yHl  set  equal  to  1,  was  then  applied  to  obtain  pre¬ 
dicted  values  of  current  density,  zpred; 


Experimental  Current  Density  (iexp)  X  103 


Fig.  2.  Parity  plot  of  ipred  vs.  iexp  parity  plot  for  the  data  in  Table  2  using  the  proposed 
Mechanism  I  Tafel-Volmer  equations  and  parameter  values  in  Table  3. 

The  resulting  values  of  ipred  were  then  calculated  for  each  polar¬ 
ization,  p exp,  in  Table  2  and  the  ipred  vs.  zeXp  results  are  plotted  in 
Fig.  2. 

For  greater  applicability  in  predicting  polarization  curves  at 
other  operating  conditions,  the  following  generalizations  are  pro¬ 
posed  to  the  above  calculation  sequence: 

(a)  The  slopes,  me ,  in  the  linear  term  of  Eqs.  (40)-(42)  are  nearly 
Arrhenius  in  their  temperature  dependence  and  can  be  approx¬ 


imately  represented  by: 

m0  »  7.3g  x  10“4  exp{2680(T)-1 )  (43) 

However,  for  greater  accuracy  at  both  interpolated  and  extrap¬ 
olated  temperatures,  the  following  are  proposed: 

me,  <303  k  *  1-475  x  10-3exp{248g(T)-1)  (44) 

mdt> 303K  «  3.78  x  1 0-4 exp{2g00(T)-1 )  (45) 

Similarly,  the  coefficients,  c#,  for  the  p335  terms  in  Eqs. 
(40)-(42)  can  be  generalized  somewhat  via  the  following: 

Ce,  <303  K  *  - 1  -36(  1 03  )T  +  4. 1 2(  1 05 )  (46) 

0?, >303 K  «  — 6.2(103)T  +  1.88(106)  (47) 

Thus,  6(p)  can  be  estimated  via  Eq.  (19),  Eqs.  (44)-(47),  and 

0  =  0o  -  mep  +  cep335  (48) 


3.3.9.  Calculation  ofipred  for  a  given  temperature  (25  °C)  and  a 
given  value  of  p 

The  calculation  sequence  proposed  in  Section  3.3.8  was  applied 
to  the  prediction  of  all  parameter  values  for  25  °C  and  an  anode 
activation  polarization,  p,  of  1  x  10-3  V.  The  applicable  equation 
numbers  and  the  results  are  summarized  in  Table  5. 

The  same  process  summarized  in  Table  5  was  applied  to  the 
prediction  of  the  polarization  curve  for  ‘Pt(l  1  0),25°C,pH2  =  1  atm, 
CM  =  1 .56  x  1 0-9  using  the  parameter  values  in  Table  5  and  applying 
Eq.  (14).  The  results  are  summarized  in  Fig.  3. 

3.4.  Analysis  of  the  published  Pt(l  10)  data  from  Protopop  off  and 
Marcus  (1985-1988) 

Protopopoff  and  Marcus,  in  a  series  of  papers,  studied  the  influ¬ 
ence  of  sulphur  on  the  electroadsorption  of  hydrogen  under  the 
H+:H2  equilibrium  potential  [20,21]  and  on  the  HER  [22]  on  a 
single-crystal  Pt(l  1  0)  surface.  Results  at  zero  sulphur  were  also 
reported  so  that  there  are  some  data  that  are  relevant  to  the 
present  paper.  Experimental  conditions  for  the  HER  study  were 
25  °C,  0.05  M  H2S04  aqueous  electrolyte,  and,  presumably,  a  pH2 
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Table  5 

Estimation  of  parameter  values  at  25  °C  and  r]  =  lx  10-3  for  Pt(l  1  0)  with  dilute 
H2S04  electrolyte  and  Cm  =  1.56  x  10-9. 


Parameter 

Source 

Estimated  value 

Note 

me 

Eq.  (44) 

6.254 

co 

Eq.  (46) 

6.72  xlO3 

UGtCH3o+,o} 

Eq.  (38) 

3.9544 

(1) 

e0 

Eq.  (19) 

0.79816 

0 

Eq.  (48) 

0.7919 

^ads 

Eq.  (31) 

9.329  xlO16 

^des 

Eq.  (32) 

4.7115  x  109 

cH2 

Eq.  (27) 

7.783  x  10“7 

[^adsCH2  (^des)  ] 

Calc’d 

3.9196 

(1) 

*'pred  (Kh2,o  —  1) 

Eq.  (14) 

8.9  x  10“5 

*'o,fwd  (Kh2,o  =  1) 

Eq.  (15) 

1.389  xl0“3 

(2) 

Io,rev 

Eq.  (16) 

1.388  xl0“3 

(2) 

(1)  For  internal  consistency,  (/<adsCH2/^des)1/2  and  {KetC^o+.ol  should  be  close  in 
value. 

(2)  For  the  parameter  values  to  show  complete  internal  consistency,  io  fwd  and  i0,rev 
should  have  identical  values. 

of  1  atm.  Current  densities  were  calculated  on  the  basis  of  the 
geometric  area  of  the  electrode.  The  Tafel  slope  for  the  zero-S  polar¬ 
ization  curve  was  0.035  ±  0.001  V dec-1  and  the  current  density  at 
0  V  was  (1.7  ±  0.05)  x  10-4Acm-2.  The  0.035  was  considered  by 
the  authors  to  be  close  enough  to  the  0.0296  value  that  would  come 
from  our  Eq.  (29)  to  support  this  HER  as  a  ‘proton  discharge  in  equi¬ 
librium  followed  by  rate-determining  chemical  recombination’  (i.e. 
our  Mechanism  I  in  reverse).  They  cited  some  additional  support¬ 
ing  literature  [23-26]  but  it  does  not  appear  that  any  of  this  work 
involved  single-crystal  Pt(l  1  0). 

As  recently  argued  [3],  a  ‘modified  Tafel  plot’  according  to  Eqs. 
(24)  and  (25)  should  be  based  on  r\  being  plotted  against  log[/(/i)— 1  ], 
not  log[i],  a  procedure  that  will  modify  the  resulting  values  of  both 
the  Tafel  slope  (i.e.  a)  and  the  Tafel  intercept  (i.e.  i0). 

The  concept  of  exchange  current  density,  the  equal  forward  and 
reverse  current  densities  at  zero  polarization  and  zero  net  current, 
implies  that  the  equilibrium  rate  of  the  HOR,  Eq.  (2),  is  the  same 
as  the  equilibrium  rate  of  the  HER,  the  reverse  of  Eq.  (2).  If  the 
mechanism,  as  argued  above,  is  the  same  for  the  HER  and  the  HOR 
in  the  equilibrium  region,  the  i0,her  and  i0,hor  should  have  the  same 
value  at  the  same  test  conditions.  Therefore  /0,hor,298k  should  also 
be  1.7  x  10-4  for  the  test  parameters  used  in  the  work  and  for  the 
geometric  area  of  the  electrode  being  used. 

The  evaluation  of  any  Tafel-Volmer  parameters  for  comparison 
with  our  results  in  Section  3.3  is  hindered  by  the  lack  of  any  mention 
of  a  value  of  CM  for  their  Pt(l  1  0)  electrode.  Since  their  i0  value 


Current  Density  (amps/cm  2) 


appears  to  be  nearly  an  order  of  magnitude  less  than  those  proposed 
in  Table  1,  either  Cm  or  some  other  experimental  parameters  were 
considerably  different  than  those  in  the  studies  by  the  Markovic 
Group. 

3.5.  Analysis  of  the  1993  published  Pt(h  kl)  data  from  Gomez  et  al. 

Gomez  et  al.  [27]  studied  the  rate  of  the  HER  in  0.5  M  H2S04 
and  “arbitrarily  took  i0  as  the  current  density  at  E  =  0”.  The  “room 
temperature”  results,  i0  (in  A  cm-2 )  and  Tafel  slope,  b  (in  V  dec-1 ), 
for  the  three  basal  Pt  orientations  were: 

Pt(l  11):  0.84  x  10-3  and  0.030 
Pt(l  0  0):  0.84  x  10-3  and  0.031 
Pt(l  1  0):  0.97  x  10-3  and  0.030 

The  almost  identical  Tafel  slope  of  ~0.03  for  all  three  crystal  ori¬ 
entations,  especially  considering  the  results  of  the  Markovic  Group 
[4-9],  probably  implies  that  the  measurements  were  concentra¬ 
tion  polarization  (i.e.  mass  transfer)  limited.  As  recently  noted 
[3],  the  concentration  polarization  expression  also  begins  with  a 
2.303RT(2F)_1  term  so  that  a  diffusion-limited  electrode  reaction 
will  show  an  apparent  Tafel  slope  of  ~9.92  x  10_5I,  from  0.0289  to 
0.0296  at  room  temperatures  from  18  to  25  °C.  For  this  reason,  the 
results  of  Gomez  et  al.  were  not  further  evaluated. 

4.  Generalizations  for  modelling  the  HOR  on  Pt(l  1 0) 

4. 1.  Modelling  in  dilute  H2S04  electrolytes 

The  implicit  assumptions  at  this  point  are  (i)  that  the  HOR  reac¬ 
tion  on  a  Pt(l  1  0)  anode  is  a  Tafel-Volmer  reaction  sequence  with 
a  ‘Mechanism  I’,  rds  Tafel  followed  by  fast  Volmer,  process;  (ii) 
that  the  dissociative  chemisorption  of  the  H2  molecule  follows  a 
Langmuir  isotherm,  and  (iii)  the  proposed  parameter  values  and 
correlations  in  Section  3.3  apply  to  the  ‘Markovic  Group’  conditions 
for  their  Pt(l  1  0)  anode:  0.05  M  H2S04,  pH2  °f  1  atm  and  CM  equal 
to  1.56x1 0-9  mol  cm-2. 

Generalizing  cH2,0  to  pH2  other  than  1  atm  requires  evaluation 
of  an  appropriate  Henry’s  Law  constant.  As  noted  in  Section  3.3.2 
with  respect  to  Eq.  (27),  Eq.  (28)  and  to  previous  work  [1  ],  given  the 
electrolyte,  the  temperature  and  the  value  ofpH2,  the  value  of  cH2,0 
can  be  estimated. 

Generalizing  to  anode  Pt(l  1  0)  catalysts  other  than  ‘Markovic 
Group  Pt(l  1  0)’  requires  adjustment  of  the  value  of  Cm-  Eqs.  (15)  and 
(16)  indicate  that  i0  varies  as  so  that  a  Pt(l  1  0)  anode  catalyst 
with  a  Cm  value  other  than  the  ‘Markovic  Group’,  ‘MG’,  value  pro¬ 
posed  in  Section  3.2  would  need  another  parameter  added,  CM, actual 
instead  of  Cm,mg»  to  properly  scale  the  predicted  i0  and  i  values.  It 
should  be  noted  that  these  CM  values  are  best  visualized  as  being 
‘per  cm2  of  geometric  Pt  electrode  area’. 

Values  of  {^etCH30+  0}  f°r  Pt(110)  when  p h2  is  other  than  1  atm 
would  also  require  adjustment.  As  noted  following  Eq.  (12),  I<et  rep¬ 
resents  the  group  {/<et,rev(ket,fwd)_1 }  where  /<et,rev  and  ket,fwd  should 
only  be  functions  of  temperature.  As  noted  around  Eqs.  (20)  and 
(21),  I<etCH30+  0  should  be  numerically  equal  to  and,  since 

JCads  should  only  be  a  function  of  temperature,  cH  0+>0  should  there¬ 
fore  vary  as  c^2,  i.e.  as  {pn^H^rV^.This  means  that,  for  a  ‘new’ 
value  of  pn2  other  than  the  Markovic  Group  “adjusted  to  a  pH2  of 
1  atm”  ‘MG’  condition, 

-1  V2 

{KetcH30+,0}pH2  ~  U4tCH3o+)0}pH2=latm[PH2HH2,MG(HH2,new)  ] 


Fig.  3.  Predicted  polarization  curve  at  25  °C  for  the  conditions  in  Table  5. 


(49) 
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Eq.  (49),  when  combined  with  the  ‘pH2  =  1  atm’  results  in  Table  3  or 
with  Eq.  (39)  for  temperatures  above  45  °C,  leads  to  the  following 
general  correlation: 

{f^etCH3Q+}pH  >45 oc  ^  2.41  x  106[pH2HH2?MG(f^H2,new)  ] 

x  exp[-3995(T)_1  ]  (50) 

The  values  of  Henry’s  constant,  Hh2  ,  in  Eq.  ( 50)  can  be  evaluated  for 
the  particular  electrolyte  as  described  elsewhere  [1  ]. 

4.2.  Generalizing  to  electrolytes  other  than  H2S04 

The  first  key  parameter  here  is  cH2  and  this  can  be  evaluated 
for  various  concentrations  of  other  acids  and  for  Nafion  membrane 
electrolytes  by  making  use  of  the  Henry’s  Law  correlations  recently 
proposed  [1  ].  As  long  as  the  pH  is  still  in  the  same  general  range  as 
that  for  the  dilute  H2S04  results,  it  can  be  assumed,  until  relevant 
experimental  data  are  available  and  appropriately  analysed,  that 
the  various  Tafel-Volmer  ‘pH2  =  1’  parameter  values  are  not  signif¬ 
icantly  changed  from  the  correlations  presented  above.  If  the  pH 
is  significantly  different  than  that  in  the  ‘Markovic  Group’  studies, 
generalization  from  the  results  in  Section  3.3  should  be  done  with 
caution. 

The  other  possibility  is  that  the  reaction  mechanism,  and  there¬ 
fore  the  values  of  the  kinetic  parameters,  differ  when  the  ionic 
makeup  of  the  electrolyte  changes.  In  this  case,  the  parameters 
derived  from  the  H2S04  studies  would  not  likely  apply. 

5.  Summary  and  conclusions 

Additional  support  has  been  provided  for  the  main  assumptions : 
(i)  that  the  HOR  reaction  on  a  Pt(l  1  0)  anode  is  a  Tafel-Volmer 
reaction  sequence  with  a  ‘Mechanism  I’,  rds  Tafel  followed  by  fast 
Volmer,  process  and  (ii)  that  the  dissociative  chemisorption  of  the 
H2  molecule  follows  a  Langmuir  isotherm  [5].  These  conclusions 
apply  to  the  ‘Markovic  Group’  conditions  for  their  Pt(l  1  0)  anode: 
0.05  M  H2S04,  Ph2  of  1  atm  and  CM  equal  to  1 .56  x  1 0-9  mol  cm-2. 

Satisfactory  Arrhenius  correlations  of  several  Tafel-Volmer 
kinetic  and  thermodynamic  parameters  for  Pt(110)  have  been 
obtained  for  the  temperature  range  from  1  °C  to  60  °C  and  for  the 
‘adjusted  to  apH2  of  1  atm’  condition  used  by  the  ‘Markovic  Group’ 
publications.  This  includes  the  parameters  {/^etC^o+l  J^ads*  ^ads*  and 
kdes.  It  is  suggested,  pending  the  availability  of  data  for  Pt(l  1  0) 
electrodes  at  temperatures  above  60  °C,  that  these  correlations  can 
be  extrapolated  for  use  at  60-90  °C,  a  typical  operating  temperature 
range  for  PEM  fuel  cells. 

Correlations  are  proposed  for  predicting  the  performance  of 
Pt(l  1  0)  electrodes  at  pH2  values  of  other  than  1  atm  and  for  CM 
values  other  than  that  recommended  by  the  Markovic  Group. 

For  modelling  purposes,  the  simplest  approach  for  current  den¬ 
sity  T  to  be  predicted  as  a  function  of  anode  activation  polarization 
p.  Alternatively,  p  can  be  predicted  as  a  function  of  current  density 
but  this  approach  would  require  an  iterative  solution.  A  calculation 
sequence  is  proposed  in  Section  3.3.8. 


There  appears  to  be  a  surprisingly  small  amount  of  published 
work  regarding  the  hydrogen  oxidation  reaction  on  Pt(  1 1  0)  anodes 
other  than  the  several  publications  of  the  ‘Markovic  Group’. 

It  is  suggested  that  these  Tafel-Volmer  parameter  values  can 
also  be  used  to  estimate  p(i)  data  for  Pt(l  1  0)  electrodes  using  other 
electrolytes  in  the  same  low-pH  region  as  the  0.05  M  H2S04  cells 
employed  by  the  Markovic  Group.  This  could  include  PEM  elec¬ 
trolytes  such  as  Nafion  although  further  work  is  needed  to  verify 
this  approach.  The  major  additional  requirement  would  be  the  pre¬ 
diction  of  Ch2  values  via  Henry’s  Law  expressions  for  the  different 
electrolytes,  a  question  that  has  recently  been  addressed  [1  ]. 

Acknowledgement 

This  work  has  been  supported  by  the  Canadian  Department  of 
National  Defence  through  funding  supplied  by  Defence  Research 
and  Development  Canada. 

References 

[1]  R.F.  Mann,  J.C.  Amphlett,  B.A.  Peppley,  C.P.  Thurgood,  J.  Power  Sources  161 
(2006) 768-774. 

[2]  R.F.  Mann,  J.C.  Amphlett,  B.A.  Peppley,  C.P.  Thurgood,  J.  Power  Sources  161 
(2006)775-781. 

[3]  R.F.  Mann,  J.C.  Amphlett,  B.A.  Peppley,  C.P.  Thurgood,  J.  Power  Sources  163 
(2007)  679-687. 

[4]  H.A.  Gasteiger,  N.M.  Markovic,  P.N.  Ross,  J.  Phys.  Chem.  99  (1995)  8290-8301. 

[5]  N.M.  Markovic,  B.N.  Grugur,  P.N.  Ross,  J.  Phys.  Chem.  101  (1997)  5405. 

[6]  N.M.  Markovic,  T.J.  Schmidt,  B.N.  Grgur,  H.A.  Gasteiger,  R.J.  Behm,  P.N.  Ross,  J. 
Phys.  Chem.  B  103  (1999)  8568-8577. 

[7]  N.M.  Markovic,  P.N.  Ross,  Surf.  Sci.  Rep.  45  (2002)  121. 

[8]  N.M.  Markovic,  The  hydrogen  electrode  reaction  and  the  electrooxidation  of  CO 
and  H2/CO  mixtures  on  well-characterized  Pt  and  Pt-bimetallic  surfaces,  2003, 
Chapter  26  of  Volume  2  of  Veilstich  et  al.  [9]. 

[9]  W.  Vielstich,  H.  Gasteiger,  A.  Lamm  (Eds.),  Handbook  of  Fuel  Cells— 
Fundamentals,  Technology  and  Applications,  vols.  1-4,  John  Wiley  &  Sons  Ltd., 
2003. 

[10]  R.  Parsons,  Trans.  Faraday  Soc.  54  (1958)  1053-1063. 

[11]  L.G.  Austen,  The  Electrochemical  Theory  of  Fuel  Cells,  in  Berger  [12]. 

[12]  C.  Berger  (Ed.),  Handbook  of  Fuel  Cell  Technology,  Prentice-Hall,  Englewood 
Cliffs,  NJ,  1968. 

[13]  J.O’M.  Bockris,  A.K.N.  Reddy,  Modern  Electrochemistry,  vols.  1&2,  Macdonald, 
Plenum  Press,  London,  New  York,  1970. 

[14]  E.  Gileadi,  E.  Kirowa-Eisner,  J.  Penciner,  Interfacial  Electrochemistry:  An  Exper¬ 
imental  Approach,  Addison-Wesley,  1975. 

[15]  T.E.  Springer,  T.  Rockward,  T.A.  Zawodzinski,  S.  Gottesfeld,  J.  Electrochem.  Soc. 
148(1)  (2001)  A11-A23. 

[16]  G.A.  Camara,  E.A.  Ticianelli,  S.  Mukerjee,  S.J.  Lee,J.  McBreen.J.  Electrochem.  Soc. 
149  (6)  (2002)  A748-A753. 

[17]  M.W.  Breiter,  Reaction  mechanisms  of  the  H2  oxidation/evolution  reaction, 
2003,  Chapter  25  of  Volume  2  of  Veilstich  et  al.  [9]. 

[18]  A.  Lasia,  Hydrogen  evolution  reaction,  2003,  Chapter  29  of  Volume  2  of  Veilstich 
et  al.  [9]. 

[19]  K.J.  Vetter,  in:  S.  Bruckensein,  B.  Howard,  Trans.  (Eds.),  Electrochemical  Kinetics, 
Academic  Press,  New  York,  1967,  pp.  516-614,  in  particular  pp.  550-551. 

[20]  P.  Marcus,  E.  Protopopoff,  Surf.  Sci.  161  (1985)  533-552. 

[21  ]  E.  Protopopoff,  P.  Marcus,  Surf.  Sci.  169  (1986)  L237-L244. 

[22]  E.  Protopopoff,  P.  Marcus,  J.  Electrochem.  Soc.  135  (1988)  3073-3075. 

[23]  J.O’M.  Bockris,  I.A.  Ammar,  A.K.M.S.  Huq,  J.  Phys.  Chem.  61  (1957)  879-886. 

[24]  S.  Schuldiner,  J.  Electrochem.  Soc.  106  (1959)  891. 

[25]  V.S.  Bagotsky,  Yu.B.  Vassiliev,  I.I.  Pyshnograeva,  Electrochim.  Acta  16  (1971) 
2141-2167. 

[26]  B.E.  Conway,  L.  Bai,  J.  Electroanal.  Chem.  198  (1986)  149-175. 

[27]  R.  Gomez,  A.  Fernandez-Vega,  J.M.  Feliu,  A.  Aldaz,  J.  Phys.  Chem.  97  (1993) 
4769-4776. 


